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SUMMARY 


Leontovich  has  proved  that  the  triangular  equilibrium 
positions  in  the  planar  Restricted  Problem  of  Three  Bodies  are 
stable  for  almost  all  admissible  mass  ratios.  It  is  shown  here 
that  the  set  of  exceptional  mass  ratios  for  which  stability 
remains  to  be  proved  or  invalidated  contains  only  one  point 


besides  the  critical  mass  ratios  of  order  two  and  three. 


1 .  INTRODUCTION 


We  consider  a  conservative  dynamical  system  with  two  degrees  of 
freedom  defined  by  the  equations  of  motion 


Zi  =  3z 


dX 


i+2 


‘i+2 


8  X 

3z , 


(i-1,2) 


with  the  Hamiltonian  function  X  analytic  at  the  origin  in  the  four 
dimensional  phase  space. 

Let  us  suppose  that  the  origin  is  a  position  of  equilibrium  for  the 
system.  Thus,  in  the  neighborhood  of  the  origin,  the  Hamiltonian  can  be 
expanded  as  a  series 


n>2 

of  powers  of  z^,  z2>  z^,  z^,  the  term  Xn  being  the  homogeneous 
component  of  the  series  with  degree  n. 

We  assume  that  there  exists  a  completely  canonical  linear  mapping 
(z1>z2,z3,z^)  (Z^.Z^Z^.Z^)  reducing  X^  to  the  form 

*2  ‘  2  “1<Z1+Z3>  +I“2(Z2+Z4) 


where  and  w2  ate  real  numbers.  Thus  the  equilibrium  is  said  to 

be  of  the  elliptic  kind,  and  the  numbers  and  to 2  represent  a  set 
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of  basic  frequencies  for  the  linear  dynamical  system  described  by  the 
Hamiltonian  function 

Given  a  fixed  integer  n  such  that  n  ^  4,  let  us  consider  the 
following  hypothesis: 

(An)  ("Restricted  Condition  of  Irrationality") — 

^1^1  ^2^2  ^  ® 

for  all  pairs  (ki*k2^  of  rational  integers  such  that 

lk! I  +  lk2l  ln* 

Under  the  condition  (A  ) ,  there  exists  at  least  one  canonical 

n 

transformation  (Z^.Z^Z^Z^)  -*  which  decomposes  the 

original  Hamiltonian  into  the  sum 


a  =  a(;1,?2,t3,c4)  -  a  +  a><n+1) 


where  3£  is  a  polynomial  in  the  variables 


T  2  2 
h  *  ?i  +  ?3- 


2  2 
*2  "  C2  +  C4 


of  the  form 


X  -  w  I  +  w2I2  +  j(AlJ  +  2BI1I2  +  CI2)  +• 


(1) 


and  where  is  a  power  series  in  C2»  ^4  beginning  with 

terms  of  degree  n  +  1  and  convergent  in  the  neighborhood  of  the  origin. 


The  coefficients  in  the  polynomial  X  do  not  depend  on  the 
integer  n  nor  on  the  manner  by  which  the  normalization  transformation 
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is  obtained.  Thus,  in  particular,  the  determinant 


D  =  det (a.  ) 


whose  elements  are  defined  as  follows 


a2  x 

^Vij/VIj=0 


dX 


ai3  "  a3i  V  31 


a33  =  °* 


(i.j-1,2), 


(i-1,2). 


is  an  invariant  of  ti  *.  Hamiltonian  X  with  respect  to  the  canonical 
transformations  leading  to  the  normal  decomposition  (1). 

At  last,  let  us  still  consider  two  more  hypotheses: 

(A  )  ("General  Condition  of  Irrationality") — For  any  pair  (k. ,k_)  of 
rational  integers,  k^w^  +  k2<AJ2  ^  0* 

(B)  The  determinant  D  is  not  zero. 

Arnol'd  (1961)  has  established  the  following  theorem.  If  the 

conditions  (A  )  and  (B)  are  verified,  then  on  each  energy  manifold 

00 

X  -  h  in  the  neighborhood  of  the  equilibrium,  there  exist  invariant 
tori  of  quasi-periodic  motions  which  divide  the  manifold.  In  consequence, 
the  equilibrium  is  stable. 

As  an  application  of  this  theorem,  Leontovich  (1962)  deduced  that 
the  triangular  equilibrium  positions  in  the  planar  Restricted  Problem  of 


-4- 


Three  Bodies  are  stable  for  all  permissible  mass  ratios  but  a  set  of 
measure  zero. 

However,  Moser  has  shown  that  Arnold's  theorem  is  still  valid 
under  the  weaker  conditions  (A^)  and  (E) . 

We  propose  here  to  apply  this  weaker  form  of  the  theorem  to  the 
triangular  positions  of  the  planar  Restricted  Problem,  in  order  to 
determine  the  mass  ratios  for  which  there  is  stability. 


2.  FIRST  ORDER 

In  the  canonical  units  and  with  respect  to  the  synodical 
barycentric  system  of  Cartesian  coordinates,  the  planar  Restricted 
Problem  of  Three  Bodies  is  described  by  the  Lagrangian  function 

£  «  y(x  +y  )  +  (xy-xy)  +  (1-p)  (—  +  j  +  y(—  +  j  P2>. 

For  the  sake  of  convenience,  we  put 

Y  -  1  -  2y. 

The  substitution  x  x  +  y/2,  y  -*■  y  +  /J/2  translates  the  origin 
of  the  coordinate  system  to  the  triangular  equilibrium  position  L^. 
There  we  expand  the  Lagrangian  function  in  power  series  of  x  and  y, 
and  we  find  (Deprit  1966)  : 

£  =  £^  +  £^  +  £^ » 
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p  1  ,.2  . 2.  ,  y  •  •  >  .  3  2  3  pr  ,9  2 

£2  =  ^(x  +y  )  +  (xy-xy)  +  jx  +  -  y/3  xy  +  j  y  , 


o  _  7  3  3  pr  2  33  2  3/^3 

*3  *  I?  **  '  16  ^  X  y  ‘  16  ^  "  16  ^  y  * 

t  37  4  25  /r  3  ,  123  2  2  ,  45  77  3  .  3  4 

4  =  ~  128  x  "I^xy+Y  xy  r  —  y/3  xy  +1237- 


Thus  to  the  first  order,  the  equations  of  motion  are  to  be  derived 
from  £ 2*  They  are 


X-2y-|-x-|  y/3y  -  0, 
y  +  2x  -  -■  y/3  x  -  I*  y  -  0. 


To  them  belongs  the  characteristic  polynomial 


4  2  27 

S  +  S  +  16  ^1_Y  *  * 


Its  roots  are  distinct  and  all  purely  imaginary  if  and  only  if 


li(l-v)  <  1/27. 


In  order  that  this  be  true,  the  mass  ratio  y  must  be  strictly  smaller 
than  Routh's  critical  mass  ratio 


1  <  f(1  "  9  ^5)  a  0.038520896... 


in  which  case,  the  four  characteristic  roots  at  the  equilibrium  are 


as  follows: 
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s^  11  lu)^  ,  82  B  lu^  i  a  — lu)^  , 

where  the  real  numbers  uj^  and  are  characterized  unambiguously 

by  the  set  of  relations 

0  <  (1)2  <  1//T  <  o)^  <  1, 

2  2 

o)i  +  «2  *  1, 

160)^2  -  27  (1-y2)  . 

These  classical  results  define  the  stability  and  the  periods  of  the 
linear  oscillations  about  the  equilibrium  omitting  the  nonlinear 

terms  in  the  Lagrangian  function.  What  interests  us  more  now  is, 
however,  to  investigate  the  perturbations  in  the  coordinates  x,y  and 
in  the  basic  frequencies  arising  from  the  presence  of  quadratic  and 
cubic  terms  on  the  right-hand  members  of  equations.  In  order  to  do  so, 
let  us  hereafter  put,  for  the  sake  of  simplicity, 

k  *  (2u)^  1/2)  “  (1  2032)  ’ 


h  = 

(9+Au)^)1/2 

(±-1,2) , 

mi  = 

1  +  4oj2 

(1-1,2), 

ni  = 

9  -  Acj2 

(1-1,2). 

Then,  following  Breakwell  and  Pringle  (1966),  we  Introduce  the  symplectic 


(2) 


matrix 
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4  = 


) 

lli.JlA 


whose  elements  are  defined  as  follows: 


all  " 

0, 

a21 

a12  = 

0, 

a22 

a13  = 

41/2kw1, 

a23 

a14  = 

-£2/2kui2 , 

a24 

a31  “ 

a41 

a32  “ 

-m2co2/2kJl2, 

a42 

a33  “ 

3y/3/2k£1o)1, 

a43 

a34  * 

-3y/3/2k£2ai2, 

a44 

-40^/k^, 

-^2/^2 ' 

-3y*/3/2kil1u)1, 

3y/3/2kJl2a)2, 

3Y/3w^/2k£,^ , 
3y/3w2/2k a2, 
n1/2kA1a>1, 
-n2/2k^2w2. 


Accordingly  the  transformation  from  the  phase  space  (x,y,p 
the  phase  space  product  of  the  angle-coordinates  (<f>^,<J>2) 
action  momenta  as  it  is  defined  by  the  equations 


x*Py^ 
and  of 


into 
the  two 


y  “  a21Qx  +  a22Q2  +  a23Pl  +  a24P2* 
px  -  a31Ql  +  a32Q2  +  a33Pl  +  a34P2’ 
py  *  a41Ql  +  a42Q2  +  a43Pl  +  a44P2’ 


wherein 
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1/2 

Q1  "  (2Ii/wi)  7  sin  *i#  (1-1,2) 

1/2 

Pi  "  (2Iic  i>  cos  *i’  (1-1*2) 


Is  completely  canonical.  Moreover,  in  the  new  phase  variables,  the 
second  order  part  of  the  Hamiltonian  assumes  the  normal  form 


Jf2  “  “  U)2I2  * 

and  its  general  solution  is 

$ .  -  oj. t  +  const. , 

( i-1,2) 

1^  =  const. 

If  the  oscillations  about  were  exactly  linear,  the  integrals  of 

motion  would,  in  fact,  be  represented  by  the  above  relations  and  the 
corresponding  orbits  would  be  given  by  the  formulae 

*1  _l/2  ,  *2  _l/2 

x  =  — -  I,  cos  <J>, - -  I„  cos 

k/2^  1  1  2  2 


y  = 


3y/3  _l/2  ,  ,  3y/3  _l/2 

— u  1  I.  cos  d>1  +  — •"  cos  4> 

k£,^/2w2  ki^^2^ 


sin  $ 


1 


4/2o)r 


k£. 


2  T1/2  ,  x 

—  I„  sin  <j> 


2' 
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3.  SECOND  ORDER  NORMALIZATION 

J.  Henrard  (1966)  has  shown  how  to  carry  on  in  a  straightforward 
manner  Birkhoff's  normalization  without  introducing  generating  functions 
and  without  inverting  power  series. 

The  coordinates  (x,y)  are  to  be  expanded  in  double  d'Alembert 
series : 


x 


(3) 


where  the  homogeneous  components  and  y^  of  degree  n  are  of  the 

form 


0<m<n 


(n-m) /2  m/2 
L1  2 


\  [C  cos (p4> ,  +q4>0 )+S  sin(p<f>..+q<j>0)] .  (4) 

n-m,m,p,q  vrYl  2  n-m,m,p,q  V^T1  My2/J 


In  (4)  the  double  summation  o  *er  the  indices  p  and  q  is  subject  to  the 
following  conventions: 

a)  p  runs  over  those  integers  in  the  interval  0  <_  p  <_  n-m 
that  have  the  same  parity  as  n-m; 

b)  q  runs  over  those  integers  in  the  interval  -m  q  <_  m 
that  have  the  same  parity  as  ra. 

In  the  developments  (3),  the  quantities  1^  and  ^  are  to  be  taken  as 
constants  of  integration,  while  and  $2  are  to  be  determined  as  linear 

functions  of  the  time  in  such  a  way  that 

*1  =  “l  +  f2n(Il,I2)’ 

n>l 

*2  "  -a)2  +  \  82n(Il,I2) 


I 
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where,  for  any  n  ^  1,  f^n  and  g2n  are  homogeneous  polynomials  of 
degree  n  in  the  actions  1^  and  I2.  The  canonical  character  of  the 
transformation  will  be  ensured  formally  by  requesting  that  the  double 
d'Alembert  series  (3)  satisfy  the  identities 


(x;y) 

■  o, 

(x;x) 

=  1, 

(y;x)  *  0, 

(x;y) 

-  0, 

(y;y)  -  1, 

(x,y)  -  0 

where  the  left-hand  members  stand  for  the  Poisson  bracket  with  respect 
to  the  phase  variables  *  thus  for  instance 

/  \  3x  9  y  3x  3y  .  3x  3y  3x  3y 

VX,y;  3^  3I1  3IX  34>1  3<t>2  3I2  3I2  3<»2  ’ 

and  where  x,y  should  be  taken  as  the  composition  derivatives 


•  i  9  x  ,  3x 

X  ”  *1  3^  *2  302  * 


* =  h  hT  +  *2  it  • 


In  this  way,  a  Birkhoff  normalizing  transformation  can  be 
constructed  entirely  by  the  method  of  undetermined  coefficients. 

As  it  is  shown  elsewhere  (Deprit  et  al  1966b),  the  homogeneous 
components  of  order  2  in  the  coordinates  x  and  y  are  solutions 
of  the  partial  differential  equations 
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7  _JL_  9  \2  3 

_lWl  3 4>1  “2  34,2  /  "  4 

(5) 


,1.0 


2(“1  3*.  ' 


(J)r 


±)  +t^lBr 


(W1  3<J>1  “2  3 4> 2  ) 


+ 


The  right-hand  members  are  the  homogeneous  components  of  order  2  obtained 
on  substituting  in  the  derivatives 


21  2 
16 


» 


_3_ 

16 


nr  Z  33 
/3x  -T 


yxy 


the  first  order  expressions  already  obtained  for  x  and  y.  From 
eliminating  in  turn  ,x  and  B^’^  from  the  system  (5),  we  arrive  at 
the  partial  differential  system 


V2B2’° 


=  <t>. 


A1A2B2 


where  the  differential  operators  on  the  left-hand  members  are 

a  /  9  9  \2  ,  2 

Arh^‘u2^)  v 

,  /  __3_  9  \ 2  .  2 

A2  "(“l  9<ji1  “2  3^  )  +  w2 


and  the  functions  on  the  right-hand  members  are  to  be  constructed  as 


/ 


follows: 
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It  is  quite  essential  to  remark  that  no  term  in  cos  <f>^,  sin  ^  or 
cos  <J>2»  sin  <^2  appears  in  X2  and  Y2 ,  so  that  $2  and  ^  are  also 
free  of  such  terms.  Hence  the  homogeneous  components  can 

obtained  in  a  straightforward  manner,  on  applying  the  differentiation 
rules : 


where 


A1A2  cosCp^  +  q<J>2>  =  Ap  cosCp^  +  qcf>2)  > 
A.^2  sin(p<^  +  q$2)  =  Ap  sinCp^  +  q4>2), 

2  2  2  2 
Ap  q  =  [  co1  -  (pu>1-qu)2 )  ]  [u>2  -  (pu^-qu^)  ]• 


Since  we  assumed  that  the  irrationality  condition  (A^)  is  verified, 
none  of  the  divisors  A2  A^  A^  Aq  2  is  zero.  Notice  that,  in 
view  of  the  inequalities  (2),  the  condition  (A^)  is  fulfilled  if 
and  only  if,  in  the  interval  0  <  p  <  the  mass  ratio  does  not  take 

the  critical  value 


M2  =  2  [  1  "  45  ]  =  °-024293897-  •  •  • 
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A  simple  sequence  of  algebraic  manipulations  leads  eventually 
to  the  components  and  ^  ^as  keen  checked  that  they 

actually  transform  the  Hamiltonian  function  into  the 

zero  function.  Both  components  are  listed  in  Table  I. 
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Table  I 


Homogeneous  components  of  order  2  in  the  coordinates  x  and 


4.  SECOND  ORDER  COEFFICIENTS  IN  THE  FREQUENCIES 


The  homogeneous  components  and  B^’^  of  the  third  order  in 

the  coordinates  x  and  y,  as  well  as  the  homogeneous  polynomials 
and  of  the  second  order  in  the  frequencies  ^  and  ^ »  satisfy  the 

partial  differential  equations 


U), 


a 

9<t>i 


-  OJ 


2  9<f), 


,1.0 


■[2(“i4:'“2it)  ^K’1 


+  2f. 


a 

a<f>i 


. .  _i_  Ri.o  _ 

1  B1 


4A} 


~  2g. 


OJ 


TT-  B^’°  +  B?*1 


2  a<j) 


•  x 


3’ 


The  right-hand  members  are  the  homogeneous  components  of  order  3 
obtained  on  substituting  in  the  derivatives 

the  expansions  obtained  for  the  coordinates  x  and  y  up  to  the  second 
order  so  far. 

As  we  did  at  the  second  order,  we  eliminate  in  turn  B^’^  and 
from  the  system  (6)  so  that  we  obtain  the  equations 


-16- 


A1A2B3*°  “  $3  "  2f2P  ~  2g2^» 


A1A2B3*1  “  V3  ~  2f2U  "  2g2V * 


where  we  put 


¥,  “ 


Q  - 


u  - 


-  f]X3  +  [2(uil  "“2  3^  +  4  y/J]Y3* 

[(ul  -  “2  -  4  ]  Y3  •  [2(“1  "  *2  -  4  Y</J]  x3’ 

_3_|\  2  9W  _3_  .1,0  .0,1,  .  ,,  _3_  .  3  yrw  3  ,0,1  ,  1,0.1 

s^i  [S  3r4}S  5*1  Bi  *Bi  )  +  (2wi^  +  r^)(“i^Bi  +  Bi  )J* 

3  2  _il  9W  _3_  „1.0  ,0,lv  .  ,  ,  3  ,  3  /r-,,  3  ,0.1  .  .1.0 

3*2  L(U12  ^2  4)(  w2  3*2  B1  B1  )  (  2w2  3*2  4  Y  3)(  w2  3*2  B1  B1 

-JL  T/  2  Jl  2w  -1  *0,1  _1,0,  _3_  3  /r-,,  3  _1,0  .0,1,1 

L(“i  "  a)(wi  Bi  +  Bi  >  "  <2“i  -  4  y/J)(“i  T^B1  "  B1  >J ’ 

2  2W  _3_  *0,1  ,  _1,0,  ,  ,  _JL  2  yrw  _3__  _1,0  _0,1 

L(w2  2  4)(_u,2  3*,  B1  +  B:  )  "  (-2w2  3*  "  4  Y  3)(~U2  3*.  B1  *  B1 

^  *  1 


_2 

3* 


We  do  not  need  to  evaluate  the  components  and  B^'^.  For 

it  will  be  sufficient  for  our  purpose  to  compute  the  coefficients  of 
cos  4>3»  sin  4^,  cos  ^  and  sin  4 2  in  the  right-hand  members  of  the 
partial  equations.  Those  are  the  critical  terms,  because  they  belong 
to  the  kernels  of  the  differential  operator  &i^2'  Now  t^iey  can  Be 
eliminated  by  a  proper  choice  of  the  coefficients  in  the  polynomials 

f 2  =  f2,0Il  +  f0,2I2  ’ 

82  —  62  o^l  «0  2*2 * 


-17- 


As  a  matter  of  fact,  it  turns  out  that  the  system  of  16  linear  equations 

in  the  four  unknowns  f ^  q»  2’  ^2  0*  2  t0  we  arr*-ved 

consistent  and  yields  a  unique  solution.  We  have  listed  it  in  Table  II. 


Table  II 

Component  of  fourth  order  in  the  normalized  Hamiltonian 


2,0 


A  =  72 


1  0^(8 1-696^+124  ojJ) 


(1-2^)  (l-5u)J) 


f 0, 2  g2,0  B 


1 

6 


2  2 

°J1°J2 


(l-2o)^)  (1-2^)  d-5^)  (1-5^) 


’0,2 


=  C  - 


x  ^(81-696^+1240)2) 

(1-20)2)  (1-5o)2) 


These  results  have  been  checked.  For  the  system  Sun-Jupiter 
(p  =  0.000953875...),  the  Birkhoff's  normalization  has  been  carried 
numerically  up  to  order  thirteen  (Deprit  et  al  1966b) ;  there  it  was 
found  that 


A  =  0.01135436, 
B  =  -0.1551412, 
C  =  1.119733. 


The  same  task  has  been  performed  (Deprit  et  al  1966a)  for  the  system 
Earth-Moon  (p  =  0.0121500...)  where  it  was  found  that 
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A  =  0.2313561, 

B  =  -1.712630, 

C  =  0.6771718. 

Both  systems  of  values  have  been  recovered  from  the  definitions  in 
Table  II  valid  for  any  permissible  mass  ratios. 

Another  check  is  provided  by  the  analytical  expansions  of  the 
two  natural  families  of  periodic  orbits  issued  from  L^,  as  they  were 
performed  by  Pedersen  (1935).  There  the  author  finds  that  the  frequency 
along  the  family  of  short  period  orbits  is  a  series 

•  2 

+  Ee  +•  •  • 

in  a  certain  orbital  parameter  e.  For  the  coefficient  E,  he  gives 
the  expression 

24  3-2007  w^-648u>^+290R(jj^-496(jj^ 

E  2  2  * 

2304w  (l-2u>p(l-5wp 


If  our  computation  is  correct,  our  coefficient  A  should  be  traced  in 
Pedersen's  coefficient  E.  But  this  is  the  case,  since 


243  -  2007w*  -  648<^  +  2908u)J  -  496io®  =  u*  (3+4wJ)  (81-696(^+124^)  ; 


hence 


E 


3+4to^ 

32w 
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Birkhoff's  normalization  cannot  be  carried  at  order  three  unless 
the  irrationality  condition  (A^)  is  assumed.  But  (/,)  means  that, 
from  the  interval  0  <  \i  <  |j^  of  permissible  mass  ratios,  besides 
the  value 

y3  =  2  (1  "  "h  =  °*013516016... 

should  be  excluded. 

Thus  A,  B,  C  are  functions  of  the  mass  ratio  on  the  interval 
0  <  m  <  out  of  which  the  critical  values  \i^  and  have  been 

taken.  The  main  qualitative  characteristics  of  these  functions  can  be 
read  from  Figure  1  in  which  we  have  plotted  respectively  A,  B,  C 
versus  the  mass  ratio. 

5.  STABILITY 

Knowing  the  fourth  order  part 

=  2(AI1  +  2BI1X2  +  ^ 

of  the  normalized  Hamiltonian,  we  compute  the  determinant  D  which 
decides  about  the  stability  at  the  equilibrium  In  the  theorems  of  Arnol'd 
and  Moser.  Thus  we  find  that 

2  2 
D  =  -(Au^  +  +  Cw^) , 

a  quantity  which  can  easily  be  expressed  as  a  rational  function  of  the 

.  „  2  2 
product  0)^2 : 
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D 


36-541tD^(j2+644oj^a)2 

d-4^2)  (4-25u)^2> 


It  is  easy  to  see  that,  in  the  interval  0  <  y  <  y^,  it  possesses  one, 
and  only  one,  zero  for  the  mass  ratio 


y  =  0.01091367...  . 
c 

Thus,  stability  of  the  equilibrium  at  cannot  be  decided  for  this 

mass  ratio  from  applying  Moser's  theorem. 

The  qualitative  features  of  the  determinant  -D  as  a  function  of 
the  mass  ratio  y  are  summarized  in  the  Figure  2. 


o 


normal 
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6.  CONCLUSIONS 

In  the  planar  Restricted  Problem  of  Three  Bodies,  the  question 
of  stability  of  the  equilateral  positions  of  equilibrium  can  be  answered 
in  the  affirmative  for  all  values  of  the  mass  ratio  y  in  the  open 
interval  0  <  y  <  y^  except  at  the  critical  mass  ratios  =  0.024293... 
and  y^  =  0.013516...  and  at  a  third  point  yc  =  0.010913...  .  At  these 
three  points,  Moser's  theorem  does  not  apply. 
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